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ABSTRACT 

Type la supernovae (SNe la) are powerful cosmological “standardizable candles” and the most 
precise distance indicators. However, a limiting factor in their use for precision cosmology rests on 
our ability to correct for the dust extinction toward them. SN 20I4J in the starburst galaxy M82, 
the closest detected SN la in three decades, provides unparalleled opportunities to study the dust 
extinction toward an SN la. In order to derive the extinction as a function of wavelength, we model 
the color excesses toward SN 20I4J, which are observationally derived over a wide wavelength range 
in terms of dust models consisting of a mixture of silicate and graphite. The resulting extinction laws 
steeply rise toward the far ultraviolet, even steeper than that of the Small Magellanic Cloud (SMC). 
We infer a visual extinction of Ay ~ 1.9 mag, a reddening ot E{B — V) ss 1.1 mag, and a total- 
to-selective extinction ratio of i?y«1.7, consistent with that previously derived from photometric, 
spectroscopic, and polarimetric observations. The size distributions of the dust in the interstellar 
medium toward SN 2014J are skewed toward substantially smaller grains than that of the Milky Way 
and the SMC. 

Subject headings: dust, extinction — galaxies: ISM — galaxies: individual (Messier 82) — supernovae: 
individual (SN 2014J) 


1. INTRODUCTION 

Type la supernovae (SNe la) are considered to be one 
of the mo st precise too ls for determining astronomical 
distances (lHowellll201lll . Because of their high luminos¬ 
ity and relatively small dispersion at the maxima of their 
bolometric light curves, they are commonly utilized as 
cosmological “standardizable candles”. The accelerated 
expansion of the Universe and the presence of dark en¬ 
ergy were disc overed through SNe la used as standard - 
izable candles ( Riess etldlllbQSHPerlmutter et al.lll999ll . 
The effectiveness of SNe la as distance indicators and 
standard candles is hampered by the systematic uncer¬ 
tainties related to their explosion mechanism and pro¬ 
genitor systems, and more importantly, the line-of-sight 
extinction. The distance d measured in parsec to a SN 
is Ig d = 0.2 {mx — Mx -|- 5 — Aa), where mx and Mx are 
its apparent and absolute magnitudes at wavelength A, 
and Ax is the extinction. As it is not easy to directly 
measure Ax-, one often measures the color excess (or red¬ 
dening) E{X — V) = Ax — Ay, where Ay is the extinction 
in the V-band (centered around 5500 A). SN reddening 
is often measured by comparing the observed SN colors 
to a zero-redden ing locus. 

ICardelli et al.l (1989; CCM) found that the Galac¬ 
tic extinction curves (or extinction laws) — the wave- 
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length dependencies of the extinction — can be closely 
parametrized by the total-to-selective extinction ratio 
Ry = Ay/E(B — V) , where the B -band centers around 
4400 A (also see iFitzoatrickl [19991 hereafter FTZ). As¬ 
tronomers often derive Ry for SNe la by fitting the ob¬ 
served E{X — V) with the i?y-based CCM formula. Once 
Ry is determined, one can apply the CCM-formula (or 
some other parameterizations) to derive Ax- However, 
we caution that the CCM- and FTZ-parameterizations 
have been derived for Galactic sightlines with 2 < Ry < 
5, and may not be valid for external galaxies. Note that 
the CCM formula is not even applicable to the Large 
and Small Magellanic Clouds (LMC, SMC: l^rdon et ^ 

iool). 


SNe la are so rare that nearby SNe la (d < 5 Mpc) are 
detected only about once a decade. SN 2014J, discov¬ 
ered in the ne arby starburst galaxy M 82 at a distance of 
d ~ 3.5 Mpc ( Dalcanton et al.l l2009ll . is the nearest SN 
la seen in the last three decades. Its proximity offers an 
unprecedented opportunity to study the extinction and 
reddening toward a SN la. The aim of this Letter is to de¬ 
rive Ry and Ax by fitting the reddening curve obtained 
by Amanullah et al. (2014) during the epoch range of 
[—5, -1-5] days around its peak brightness (iJ2]) using the 
silicate-graphite model (SJ3]). The results are presented in 
§3, discussed in §4, and summarized in §5. 


2. COLOR-EXCESS CURVES OF SN 2014J 


Various studies have been carried out to determine the 
Ry value for the sigh t line toward SN 2 0 14J (e.g., see 


Amanullah et al. 

2014: 

Folev et al. 

2014 

1 " 

Goobar et al.l 

2014 

Marion et al. 20151 Weltv et al.l 

2014). More 


specifically, based on the UV to near-IR photometry of 
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Fig. 1.— (a) Comparison of the color-excess curve of F14 (red 
solid line) with that of A14 (blue triangles). The blue dashed line 
is the best-fit FTZ model of A14, while the green dotted line is 
the best-fit two-component model (i.e., a combination of interstel¬ 
lar reddening and circumstellar scattering) for all epochs (“CSMD 
All”) of SN 2014J of F14. (b) Comparison of the color-excess curves 
of F14 (red solid line) and A14 (blue triangles) with our model pre¬ 
dictions (magenta dashed line for F14, and cyan dot-dashed line for 
A14). (c) Residuals between our model color-excess E(V — A)™°^ 
for A14 and (i) the data points of A14 (blue triangles), (ii) the 
FTZ fit of A14 (blue dashed line), (iii) the color-excess curve of 
F14 (red solid line), (iv) the CSMD fit of F14 (green dotted line), 
and (v) our model color-excess for F14 (magenta dashed line). 


SN 2014J obtained with the WFC3 filters aboard the 
Hubble Space Telescope (HST) and ground-based tele¬ 
scopes, lAmanullah et al.1 (2014, hereafter A14) deter¬ 
mined the reddening curve E{X — V) for SN 2014J over 
a wide wavelength range by comparing the colors of 
SN 2014J with that of SN 2011fe, a reddening-free SN la. 
They derived Ry ~ 1.4 zb 0.1 by fitting the observation- 
ally determined E{X — V) with three template extinction 
laws: an MW-type law as parameterized by CCM or FTZ 
(see FigureUb), a th ird-order polynomial SALT2 law (see 
iBetoule et ^.1 20141) , and a power-law parameterization 
~ which was shown to be a Rood appro ximation 
for multiple scattering scenarios ( Goobaijl2008li . 

We note that A14 shifted the effective wavelengths of 
the HST/WFC3 filters (see their Figure 3). In the follow¬ 
ing, we adopt the wavelength-shifted color-excess data of 


A14. 

Foley et ah (2014, hereafter F14) determined the ex¬ 
tinction curve toward SN 2014J at t = —6.4, -1-24.1 d by 
comparing the UV/near-IR SED of SN 2014J with that 
of SN 201 Ife. As shown in Figure [T] the color-excess 
curve of F14 closely resembles that of A14. The differ¬ 
ence between the F14 curve and that of A14, on aver¬ 
age, is ^0.19mag (see Figure[TJ:). Using FTZ and CCM 
models, F14 derived Ry ~ 1.66 ±0.03 and r; 1.44 ±0.03 
for SN 2014J, respectively. F14 also argued that a two- 
component model (CSMD) with both a circumstellar 
scattering component of A\ ^ and an FTZ in¬ 

terstellar reddening component of Ry ~ 2.6 could best 
account for the observed properties of SN 2014J. 


3. DUST MODEL 


We consider the silicate-graphite grain model that con¬ 
sists of amorphous silicate and graphite ( Draine fc Led 
Il984!) . We adopt the same exponentially cutoff power- 
law size distribution for both components: drii/da = 
e:xjp {—a/ab) for the size range of 50A < a < 
1 /rm, where a is the spherical radius of the dust, nn is 
the number density of H nuclei, drii is the number den¬ 
sity of dust of type i with radii in the interval [a, a ± da], 
a and Of, are, respectively, the power index and expo¬ 
nential cutoff size, and Bi is a constant related to the 
total amount of dust of type i. The total extinction at 
wavelength A is given by 


Ax/Nh = (2.51oge)y^ 


1 dn • 

da—^C'ext,Aa,A) , (1) 

riH da 


where the summation is over the two grain types. Ah 
(tih) is the H column (number) density, and ( 76 x 1 . 1 ( 0 , A) 
is the extinction cross section of grain type i of size a 
at wavelength A ca lculated from Mie t heory using the 
optical constants of iDraine fc Led (Il984l) . 

For a given set of parameters a and ab, we derive the 
constant Bi from the abundances of the dust-forming ele¬ 
ments. Let [X/H]jg,^^ be the total interstellar abundance 
of element X (i.e., Fe, Mg, Si, O, and C) relative to H in 
the interstellar medium (ISM) of M82, [X/Hjg,^^ be the 
amount of X in the gas phase, [X/H]jj^j,^ be the amount of 
X contained in dust (obviously, = [X/H] 

ISM 

[X/Hjgj^g), and fix be the atomic weight of X. Let psii ~ 
3.5gcm“^ and pgra ~ 2.24 g cm”^ , respectively, be the 
mass density of amorphous silicate and graphite. For a 
chosen set of dust depletions, we derive Bi from the dust 
size distributions: 

PFe [Fe/H]j^^t ± /iMg [Mg/H]j„,^t ± psi [Si/H]^^^^ 


nnBsii = 


f da (47r/3) a^ Psn a “ exp (—a/ab) 

4 X po 




Mc [C/H] 


<iust 


( 2 ) 

(3) 


/ da (47r/3) Pgra a “exp(—a/ob) 
where we assume a stoichiometric composition of 
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Mg 2 xFe 2 (i_x)Si 04 for amorphous silicate. 

M82 is a prototypical starburst galaxy, experiencing 
a major star formation episode in its n uclear region, 
with s trong superwind and SN activity. lOriglia et al.l 
( 2004ll obtained the stellar abundances in the nuclear 
region of M82, and compared them with those of the 
hot gas derived from the nuclear X- ray spectra. Com¬ 
pared with the solar abundance of iGrevesse fc SauvaJ 
( 1998ll . both the cool stars and the hot gas in M82 sug¬ 
gest a reduction of Fe/H by ~ —0.35 ± 0.2 dex (i.e., 
[Fe/H]jg]yj 14.1 ppm) and an overall Si/Fe and Mg/Fe 
enhancement by ~0.4 and 0.5 dex, respectively (i.e., 
[Si/H]jg^ ~ 35.5 ppm, [Mg/H]jgj^ r; 44.7 ppm). Oxy¬ 
gen is enhanced by ^0.3 dex in stars and reduced by 
^^0.2 dex in the hot gasB T he stellar abundance of C 
derived bv lOriglia et al.l (l2004h is only ^ 1/4 of solar (i.e., 
[C/H]ism R:! 83.2 ppm). 

Similar to the Galactic ISM, we assume in M82 that 
Fe, Mg and Si are all locked up in silicate dust (i.e., 
[Fe/H]^,,, [Fe/H]igj^ 14.1 ppm, [Mg/H]^,,, r. 
[Mg/H]igM - 44.7pp m^/H]j'^^, « [Si/H]igM ~ 

35.5 ppm; lOriglia et al.l l2004ll . For carbon, it is less 
clear. In the Galactic ISM, a substantial fraction (« 
42%)0 of the total car bon abundance is in the gas phase 
([G/H]g,,^j, Ri 140 ppm. iGardelli et al.lflQ^ . For carbon 
dust (i.e., graphite) in M82, we will consider three cases 
[G/H]^ = [C/Hlyg^ Ri 83.2 ppm (i.e., all G is locked 

up in dust), = 1/2 [G/H]jg^j Ri 41.6 ppm, 

and [G/H]^,jg^ = Oppm (i.e., all G is in the gas phase 
and the dust model only consists of amorphous silicate). 
With [G/H]j^gj = 41.6 ppm, the model has a silicate-to- 
graphite mass ratio of Wsii/mgra ~ 10. This is close to 
that of the SMG (wsii/mgra ~ 12, Li et al. 2006). 

To facilitate a direct comparison with the color ex¬ 
cesses E{y — A)°’^® derived by A14 for SN 2014J, we first 
calculate A\ from EqlT] and then convert to reddening 
i5(G —A)™°'^ = Av — A\. For simplicity, the model color- 
excess has not been convolved with the HST/WFG3 fil¬ 
ters. We evaluate the goodness of fitting by 

x2 [E{V - A,-)”^°^ - E{V - 

d.O.f -^obs -^para 

(4) 

where E{V — is the observed color excess toward 

SN 2014J at wavelength Xj derived by A14, cr(Aj) is the 
uncertainty of E{V — E{V — Aj)^°^ is the model 

color excess at Aj, A^obs = 16 is the number of observa¬ 
tional data points, and A^para is the number of adjustable 
parameters. 


4. RESULTS 


^ An accurate knowledge of the 0/H abundance is not re¬ 
quired. The amount of O/H locked up in dust is controlled 
by Si/H: = 4[Si/H]^^g^ for a silicate composition of 

Mg2xFe2(i_x)Si04. 

^ We take the Galactic interste llar carbon abundance to be solar: 
[C/H]jgjyf [C/H]q 331 ppm jGrevesse fc Sau\^ll998l^ . 



Fig. 2.— Fitting the color excesses between —5 and -|-5 days 
from the B maximum of SN 2014J (blue triangles; A14) with the 
silicate-graphite model (red solid lines), assuming all (a), half (b), 
and none (c) of the carbon elements are depleted in graphite. The 
dotted lines plot the GCM reddening curves with the Ry values 
corresponding to that of the models (a: Ry = 1.7, b: Ry = 1.S, 
c: Ry = 1.7). For comparison, the Galactic average of Ry =3.1 
(dashed line) and CGM reddening curve of Ry = 1.4 (dot-dashed 
line) are also shown. 


In fitting the color excesses of SN 2014J observation- 
ally determined by A14, we have three parameters (i.e., 
A^para = 3): Of, Ub, End Vh. As shown in Figure^ excel¬ 
lent fits to the observed color excesses can be achieved 
by varying a and ab for a given [G/Hj^^g,. (see Tabled]). 

The best-fit models derive Ay Rs 1.9 mag, E{B — V) ^ 
1.1 mag, and Ry R:! 1.7. The reddening and Ay are con- 


1.3 mag and Ay ~ 1.8-2.0 mas 

r (lAmanuilah et al.l 

2014 

Ashall et al. 

2014 iFolev et al. 

2014: Goobar et al.l 

2014 

Weltv et al.l 

20141). while the model Ry values are some- 


what larger than most of the earlier results , i.e., Ry ^1.4 


(lAmanuilah et al.ll20l4 feown et ai]l2015l : [Goobar et al.l 


12014 : Marion et al.l 2015 ). 

The best-fit models suggest TVh ~ 1.3-3.3x 10^^ cm“^ 
for the sightline toward SN 2014J, somewhat higher than 
that of the hot gas in M82, Vh ~ (7. 9 ± 0.7) x 10^| cm"^ , 
derived from Chandra observations (jOriglia et al.l[2003l . 

In Figure[31 we show the best-ht models shown in Fig¬ 
ured) in terms of Ax/Ay. For comparison, we also show 
the GGM reddening curve of Ry = l.A, the average ex¬ 
tinction curves for the MW {Ry = 3.1) and the SMC 























































4 


TABLE 1 

Model Parameters and Results 


(ppm) 

'^H^gra 


Q 

(lb 

ifim) 

ATh 

(1022 cm-2) 

x2/d.o.f 

£;(B - V) 

(mag) 

Av 

(mag) 

Rv 

83.2 

1.3 X 10“*® 

4.0 X 10-2® 

2.5 

0.03 

1.3 

0.45 

1.1 

1.9 

1.7 

41.6 

1.0 X 10-2® 

1.6 X 10-24 

3.2 

0.05 

1.9 

0.47 

1.1 

1.9 

1.8 

0 

0 

1.2 X 10-2® 

3.7 

0.09 

3.3 

0.53 

1.1 

1.9 

1.7 

83.2*> 

3.7 X 10-23 

1.4 X 10-23 

3.1 

0.05 

1.6 

0.51 

1.1 

2.0 

1.8 

41.6*' 

2.5 X 10-24 

4.7 X 10-23 

3.3 

0.06 

2.1 

0.51 

1.1 

1.9 

1.8 

83.2'= 

4.7 X 10-*® 

1.4 X 10-1® 

2.4 

0.03 

1.3 

0.43 

1.1 

1.9 

1.7 

83.2'* 

3.0 X 10-21 

9.1 X 10-22 

2.8 

0.03 

1.5 

2.96 

1.1 

1.9 

1.7 


^Assumed carbon depletion in graphite or amorphous carbon. 
^Amorphous carbon. 

‘^The F218W and F225W data points of A14 were excluded. 
■^Best fit for F14 at t = —6.4d. 
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Fig. 3. — Comparison of the model extinction curves of SN 2014J (red solid) with that for the CCM Ry = 3.1 (black dash-dotted), Ry = 1.4 
(green dotted), Ry = 1.7 or 1.8 (blue dashed), the SMC bar (magenta short dashed), and the Calzetti attenuation law for starbursts (cyan 
dash-dot-dotted). 


bar. The model extinction laws for SN 2014J all exhibit 
a rapid far-UV rise which is even much steeper than that 
of the SMC bar. Unlike the SMC bar, the SN 2014J 
extinction laws display an appreciable extinction bump 
at 2175A for > 0. We also show in Figure^ 

the extinction curves predicted from the CCM formula 
with the corresponding model-derived Ry values (i.e., 
1 .8). It is seen that they substantially differ 
from that calculated from the dust models. 

As mentioned in SJ51 A14 shifted the effective wave¬ 


lengths of the HST/WFC3 filters, especially for the 
bluest F218W and F225W bands, which are highly de¬ 
pendent on the SN spectrum and the reddening law. To 
examine the effects of the wavelength-shifts, we have also 
modeled the observed color-excess curve of A14 by ex¬ 
cluding the F218W and F225W data points (which could 
be most affected). The results do not show any appre¬ 
ciable differences (see Table [T]). We have also modeled 
the F14 curve (at t = —6.4d). As shown in Figure [TJj 
and Table [U the model extinction curves and the result- 
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a ((xm) a (ixm) a (|xm) 


Fig. 4.— Comparison of the size distributions of silicate (blue solid) and graphite (red solid) derived for SN 2014J with that of the 
MW Rv = 3.1 and SMC (silicate: cyan dotted; graphite: magenta dashed; see WDOl). Also shown is the MRN size distribution (black 
dot-dashed). 


ing Av and Rv are very close to that derived from the 
A14 curve. Figure [U: shows the differences between the 
best-ht model color-excesses for A14 and that of the FTZ 
model for A14, the CSMD model for F14, as well as the 
best-ht model for F14. They are generally within the ob¬ 
servational uncertainties. However, as shown in Figure^ 
the model extinction at A < 0.3 /im differs considerably 
from the CCM parameterization. The latter is known to 
be invalid for extragalactic sightlines. 

In Figure 0] we show the dust size distributions 
derived from the models with = 83.2, 41.6 

and Oppm. Co mpared with that of the MW aver- 
age of Ry = 3.1 iMathis et ahlflQT^ (hereafter MRN); 
IWeingartner fc Draind 200ll (hereafter WDOl)], the size 
distributions of the dust in the ISM toward SN 2014J 
are skewed toward substantially smaller grains. The 
MRN size distribution is a power-law dn/da oc a~^'^ 
in the size range of 50 A < a < 0.25 /rm for both dust 
components. The WDOl size distributions extend the 
lower cutoff size to amin = 3.5 A with the smallest grains 
(a 5-, 50 A) constrained by the near- and mid-IR emission 
fsee lLi fc Drainel [20011) . The dust model presented here 
for SN 2014J assumes Omin = 50A since the UV extinc¬ 
tion ca nnot constrain the exact size of nano-sized dust 
('see lril [20041) . 


5. DISCUSSION 


5.1. The Extinction Curves 


The host galaxy of SN 2014J, M82, is regarded as 
one of the archetypical s tarburst and superwind galax¬ 
ies. ICalzetti et al.l ( 19941) derived the internal dust ex¬ 
tinction in starbursts from their overall emission spec¬ 
tra. The inferred attenuation curve is much flatter than 
that derived for SN 2014J (see Figure [3]). While the SN 
2014J model extinction curves for > 0 exhibit 


an appreciable bump at 2175 A (e.g., AT 2175 « 1.8 for 
the [C/H]^^g^ Ri 83.2 ppm model), the starburst attenu¬ 
ation curve shows no evidence for the 2175 A bump. It 
is not clear to what extent the flatness of the apparent 
starburst attenuation curve may be due to the effects of 
radiative transfer in optically thick distributions of stars 
a nd dust. _ 

iHutton et al.l (l2014h analyzed the UV images of M82 
taken by the UV/Optical Telescope (UVOT) on board 
Swift. The color-color diagram obtained with the UVW2 


(2033 A), UVM2 (2229 A), and UVWl (2591 A) filters 
is especially sensitive to the presence of the 2175 A 
bump. They examined the color-color diagram and ar¬ 
gued against a “bump-less” Calzetti-type law. 

In the silicate-graphite model presented here, 
the 2175 A extinction bump is produced by small 
graphite dust. I f we consider amorphous carbon 
( Rouleau fc Martini Il99lh instead of graphite, the 
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2175 A bump will be absent. Excellent fits can also 
be obtained from a mixture of silicate and amorphous 
carbon (see Table [Ij. Unfortunately, neither the A14 
color-excess data points nor the F14 extinction curve 
covered the 2175 A bump. Due to the lack of spectral 
features (e.g., the 2175 A bump), we are not able to 
constrain the exact composition and quantity of the 
carbon dust component (see Table [T]). 

The extinction laws derived for SN 2014J are even 
steeper than that of the SMC bar (see Figure [3]); corre¬ 
spondingly, the dust sizes of SN 2014J are smaller than 
that of the SMC bar (see Figure S]). This may be related 
to the intense UV radiation and shocks associated with 
star formation in M82 that could destroy the dust and 
lead to a predominance of small grains. 

In addition to the 2175 A bump, the models presented 
here also predict two absorption features around 9.7 and 
18 ^m arising from amorphous silicate (see FigurelS]). De¬ 
pending on [C/Hjjjjg,., the optical depth of t he 9.7 Mm fea¬ 
ture ( Ath.?) ranges from ^0.34 to ~0.86. iBeirao et al.l 
(|2008ll reported the detection of the 9.7 /rm feature in the 
Spitzer/IKS spectra of the central region of M82 with 
an optical depth of Arg.r 0.3-3.1. iTelesco et al.l ( 2015ll 
obtained the 8-13 /rm mid-IR spectra of SN 2014J at 
57 to 137 days after explosion and did not detect the 
9.7 p,m feature. 


5.2. Rv 

It is often suggested that the extinction laws toward 
SNe la are “non-standard” or “unusual” in the sense that 
unlike the MW mean value of Rv ~ 3.1, SNe la often 
have a much smaller Ry (see Figure 3 of Howell 2011), 
indicating steep UV extinction. 

There are several examples of highly reddened SNe la 
for which Ry c an be measured direct l y. They all have 
Ry < 2 (e . g., Krisciunas et al.l 20061: Elias-Rosa et al.l 
20061 120081: iNobili fc Goobad 120081: IWang et al.l 12008 


Folatelli et ai.r 2010ll . However, many low reddened 


[E{B — V) < 0 .3 mag] ones have Ry values clo se to 
that of the MW ( Mandel et al.l[20TTI : |PhilliDsll2012ll . On 


SNe la have found Ry < 2 (e.g.. Nobili et al.l 

2005 

Guv et al. 

2005: Hicken et al.l 2009 

; [Folatelli et al.l 

2010 

Burns et al. 201 ill ® 

Betoule et al. 

(20141 analyzed 740 


low- and high-z SNe la and yields Ry ^ 2. 

Many lines of evidence show that the reddening law 
to SN 2014J has a l ow value of Ry » 1.4 (see Q. 
iPatat et"^ ( 2015ll and iKawabata et al.l (|2014ll presented 


spectropolarimetric and optical/near-IR multi-band po- 
larimetric observations of SN 2014J and both indicated 
a low value of Ry < 2. However, relatively larger Ry 
values have also been suggested for SN 2014J. Using 
the equivalent wid ths of 10 diffuse interstellar bands, 
IWeltv et al.l ( 2014ll yielded E{B — U) ^ 0.71 ± 0.11 mag 


® However, it has also been suggested that the low values of 
Ry derived from large samples may partly result from the poor 
ass umptions about the intrinsic color dis tribution of SNe la (e.g., 
see iFolev et al.ll201ll : iMandel et al.ll201in . 


and Rv ~ 2.7 for SN 2014J. 

The extinction curves Ax/Ay and Ry values derived 
here for SN 2014J are based on detailed dust modeling 
of the observed color excesses. They are generally con¬ 
sistent with previous studies of Ry < 2 for SN 2014J. 
However, we caution the use of Ry to derive an extinc¬ 
tion law for SN 2014J (or any extragalactic sightlines) 
since, as demonstrated in Figure [31 the model extinction 
curves differ substantially from that calculated from the 
CCM formula. 


5.3. Interstellar or Circumstellar? 

It has been suggested that multiple scattering by 
circumstellar dust surrounding their progenitors could 
explain the non-standard reddening o b served in the 
lines of sight to SNe la ( Wand l2005t iGoobail 120081 : 
lAmanullah fc Goobaill201lh . However, F14 argued that 
the wide range of observed properties for SN 2014J may 
be caused by a combination of interstell ar reddening and 
scatte ring off circumstellar material. [Johansson et al.l 
(201^ analyzed the 3.6 and 4.5/rm Spitzer/IKAC data 
of SN 2014J and detected no significant IR excess. They 
hence placed an upper limit of Mdust ^ 10“® Mq on the 
pre-existing dust in the circumstellar environment of SN 
2014J. This is insufficient to accou nt for the observed 
non-standard reddening. Moreover. [Brown et al.l ( 2015[ l 
analyzed the light curves and color evolution obtained 
with Swift/\]WOT. They argued that these observations 
are inconsistent with a contribution scattered into the 
line of sight by circumstellar dust. 


6. CONCLUSIONS 

The extinction toward SN 2014J in M82 is derived 
as a function of wavelength from fitting the observed 
color excesses with a mixture of silicate and graphite or 
amorphous carbon dust. Insensitive to the exact car¬ 
bon dust composition and quantity, the model derives 
Ay Ri 1.9 mag, E{B — V) ~ 1.1 mag, and Ry ~ 1.7, gen¬ 
erally consistent with those reported in the literature. 
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11173007, 11373015, 11178003, and 11325313, 973 Pro¬ 
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1109039, NNX13AE63G, and the Fundamental Research 
Funds for the Gentral Universities. 
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TABLE 2 

Modeled Extinction towards SN 2014J® 


= 41.6ppm = 81.3ppm 


Band 

A 

(/im) 

Aa 

(mag) 

C*ext 

(cm^/H) 

(mag) 

C*ext 

(cm^/H) 

Ly edge 

0.091 

39.87 

2.06x10-21 

29.32 

2.22x10-21 

Ly« 

0.122 

29.40 

1.52x10-21 

22.85 

1.73x10-21 

UVW2/UVOT 

0.203 

13.87 

7.16x10-22 

13.09 

9.92x10-22 

UVM2/UVOT 

0.223 

12.53 

6.47x10-22 

12.44 

9.43x10-22 

UVWl/UVOT 

0.259 

8.64 

4.46x10-22 

8.63 

6.54x10-22 

F225W/HST 

0.287 

6.98 

3.60x10-22 

6.96 

5.28x10-22 

F275W/HST 

0.290 

6.84 

3.53x10-22 

6.82 

5.17x10-22 

F218W/HST 

0.311 

5.97 

3.08x10-22 

5.96 

4.52x10-22 

F336W/HST 

0.340 

5.03 

2.59x10-22 

5.02 

3.81x10-22 

«/SDSS 

0.355 

4.63 

2.39x10-22 

4.62 

3.50x10-22 

U 

0.365 

4.39 

2.27x10-22 

4.39 

3.33x10-22 

F438W/HST 

0.433 

3.15 

1.63x10-22 

3.16 

2.40x10-22 

B 

0.440 

3.06 

1.58x10-22 

3.06 

2.32x10-22 

F467M/HST 

0.468 

2.70 

1.39x10-22 

2.71 

2.05x10-22 

3/SDSS 

0.469 

2.69 

1.39x10-22 

2.70 

2.05x10-22 

V 

0.550 

1.94 

1.00x10-22 

1.94 

1.47x10-22 

F555W/HST 

0.550 

1.94 

1.00x10-22 

1.94 

1.47x10-22 

r/SDSS 

0.617 

1.54 

7.93x10-23 

1.53 

1.16x10-22 

F631N/HST 

0.630 

1.48 

7.62x10-23 

1.47 

1.12x10-22 

R 

0.700 

1.19 

6.13x10-23 

1.18 

8.98x10-23 

j/SDSS 

0.748 

1.04 

5.35x10-23 

1.03 

7.82x10-23 

F814W/HST 

0.792 

0.92 

4.75x10-23 

0.92 

6.95x10-23 

F845M/HST 

0.863 

0.77 

3.99x10-23 

0.77 

5.83x10-23 

z/SDSS 

0.893 

0.72 

3.72x10-23 

0.72 

5.44x10-23 

I 

0.900 

0.71 

3.66x10-23 

0.71 

5.35x10-23 

J/2MASS 

1.235 

0.38 

1.94x10-23 

0.37 

2.82x10-23 

H/2MASS 

1.662 

0.22 

1.14x10-23 

0.22 

1.64x10-23 

Ks/2MASS 

2.159 

0.15 

7.53x10-24 

0.14 

1.07x10-23 

Wl/WISE 

3.353 

0.08 

4.32x10-24 

0.08 

5.83x10-24 

L 

3.450 

0.08 

4.18x10-24 

0.07 

5.63x10-24 

[3.6]/IRAC 

3.545 

0.08 

4.06x10-24 

0.07 

5.45x10-24 

[4.5]/IRAC 

4.442 

0.06 

3.18x10-24 

0.05 

4.13x10-24 

W2/WISE 

4.603 

0.06 

3.07x10-24 

0.05 

3.97x10-24 

M 

4.800 

0.06 

2.96x10-24 

0.05 

3.81x10-24 

[5.8]/IRAC 

5.675 

0.05 

2.78x10-24 

0.05 

3.42x10-24 

[8.0]/IRAC 

7.760 

0.10 

5.20x10-24 

0.07 

5.60x10-24 

N 

10.600 

0.40 

2.06x10-23 

0.28 

2.08x10-23 

W3/WISE 

11.561 

0.27 

1.40x10-23 

0.19 

1.42x10-23 

Q 

21.000 

0.16 

8.10x10-24 

0.11 

8.34x10-24 

W4/WISE 

22.088 

0.14 

7.21x10-24 

0.10 

7.47x10-24 


^The extinction results with = 41.6ppm are recommended. Because of the space limitation of ApJL, this table is not shown in 

the published version of this Letter. 
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